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Due to nontrivial orbital pairing symmetry, surfaces and interfaces of high-T c superconductors 
support states which violate time-reversal (T-) symmetry. Such naturally degenerate states, useful 
as working states of a qubit, are standard for atomic or molecular-size qubit prototypes (e.g. based 
on nuclear spins), but exceptional for mesoscopic qubits. (In particular, they hold promise of a 
better scalability.) In these lectures I review the physics of T-breaking on surfaces and interfaces of 
high-T c superconductors; then describe existing proposals for high-T c based qubits and the current 
state of experiments; finally, I discuss the decoherence sources in the system, open questions, and 
future research directions. 



I. TIME-REVERSAL SYMMETRY BREAKING ON SURFACES AND INTERFACES OF HIGH-T C 

SUPERCONDUCTORS 

A. Description of transport in high-T c structures. 

The lack of accepted microscopic theory of superconductivity in high-T c cuprates did not prevent successful research 
in this field since 1986. We have in several respects the repetition of the situation ca 1938, but with a clear advantage 
of already having BCS theory to provide insight and language for phenomenological treatment. 

On this level, high-T c superconductors are successfully described by Gor'kov equations for normal and anomalous 
Green's functions 1 , which in Matsubara representation are defined in the usual way: 

G aP (k, r; k', t') = - (T T a a {T)al(T')) , 

F af3 = (k, r; k', r') = (T T a kQ (r)aL k ^(r')) , (1) 

F+ fj = (k,r;k',r') = (r T a^ a {r)a{ , p {r')) . 

The only difference from conventional superconductivity is in the nontrivial symmetry of the pairing potential, 

Hint =2 ^ Va ' 3 ' A ^( k ' k ') a -k+q/2. Ct a k+q/2,/3 a k'+q/2,Aa-k'+q/2, A1 - ( 2 ) 
k,k',q 

As usual, we write equations of motion for each of these functions and use the "anomalous mean field" recipe to 
decouple the four-operator products, (aWaa) — > (aW) (aa). This brings out modified self-consistency relations for 
the order parameter, 

A Q/3 (k, q) = - ]T Vp aM (k, k')F v (k' + q/2, r; k' - q/2, r), 

k' 

A+ i (k, q )=-^^, M (k',k)F ( +(k'-q/2,r;k' + q/2,r), (3) 

k' 

and Gor'kov equations, which in the spatially uniform, stationary system read: 

(iuj n - Cfe)G ! ai a(k,a;„) + A Q7 F^(k, u n ) = 5 a p\ 

(tw n + &)i ; * 9 (k,w n ) + A+ 7 G 7/3 (k,^„) = 0; (4) 
(iuj n - £ fc )F Q/3 (k,w„) - A Q7 G / 3 7 (-k, -u n ) = 0. 



* Formerly also at: D-Wave Systems Inc., 320-1985 W Broadway, Vancouver, B.C., V6J 4Y3, Canada. 
^Lectures at the summer school "Quantum Computation at the Atomic Scale", Istanbul, 2003; Turk. J. Phys. 27 (2003) 491; references 
updated. 
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Here is the fourier transform of the kinetic energy operator, £ = (2m) _1 p 2 — //, /i being the chemical potential. 

It was established, that in high-T c cuprates, like YBCO, the order parameter is a spin singlet with d-wavc orbital 
symmetry, 

A Q/3 (k) = 6 a pA(k), A(k) cx cos 2 (k x ) - cos 2 (k y ) (5) 

(with axes chosen along crystallographic directions (1,0,0), (0,1,0) in the cuprate layer )Ai*i2ii 3 .. For the following, the 
most important consequence of this symmetry is the sign change of A(k) for certain directions. This means, first, that 
there exists an intrinsic phase shift of 7r between different directions in the crystal; second, that in certain (nodal) 
directions the order parameter is zero, and therefore the quasiparticle excitation spectrum is not gapped. 

On the spatial scale exceeding the coherence length, £o, it is more convenient to use the Eilenberger equations 2 , 
which follow from Q in the quasiclassical limit. This is certainly justified in high-T c cuprates with their small £o- 

The Eilenberger equations are conveniently written in matrix form, 

v f • GK) + [w„ts + A, GK)] = 0. (6) 
Here the matrix Green's function and order parameter, 

G(v F ,r;c„) = ( ) ; A(v F ,r;c„) = ( ° , A ( v *' r ) ) , (7) 

depend both on position r and on (direction of) Fermi velocity v F . (The layered structure of high-T c cuprates allows 
us to reduce the problem to two dimensions to a good accuracy, unless we have to consider, for example, a twist 
junction 3 , or tunneling in the c-direction.) The components of G, obtained from Gor'kov's functions by integration 



over energies, satisfy the normalization condition, g Un = y 1 — fu n fu„> an d the self-consistency relation becomes 

A(v F , r) = 2^iV(0)T £ (V(v F , v' F )f Un (v F , r)) g . (8) 

u„>0 

Here the angle averaging ()g = f Q * 

In a little different language, the same results are obtained with the Andreev approximation in the Bogoliubov-de 
Gennes equations for the components of the single-bogolon wave function, (u(r), v(r)) T . The original Bogoliubov-de 
Gennes equations are obtained in the process of diagonalization of pairing BCS Hamiltonian*: 

^V 2 - n A fe (r) ^ ( u k (r) \ ^ ( -y F ■ V A fe (r) \ f u k (r) \ _ „J u k (r) 



Here the quasimomentum k labels the bogolon state, E k is the excitation energy, and the self-consistency relation 
reads 

A£(r, T) = V(k, k')u* k , (r)v k , (r) tanh Ek ' ^ (r)] . (10) 

k' 

Unlike the Gor'kov equations Q (or the initial Bogoliubov-de Gennes equations), the equations (l(ilf)l) are of the 
first order in gradients, which allows us to introduce quasiclassical trajectories (characteristics) along v F and solve 
the corresponding equations by integration along these trajectories, with proper boundary conditions. 

It is known that a quasiparticle (electron or hole, described, in terms of Bogoliubov-de Gennes equations, by a 
vector ipk = (1, 0) T exp(ikr) ((0, 1) T exp(zkr)) respectively) impinging on the superconductor from the normal metal 
can undergo an Andreev reflection, switching the branch of the excitation spectrum, acquiring an additional phase, 
and almost exactly reversing the direction of its group velocity (this happens because for an electron and a hole with 
the same momentum k group velocities are opposite): 



where 



tp^1l A -^, (11) 



n A = \ „-m/2-i X n ) ■ (12) 



The phase x is the phase of the superconducting order parameter; the (— 7r/2)-shift is exact in the limit when the 
quasiparticle energy is much less than the superconducting gap (generally it is some energy-dependent function 5(E)). 
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Now consider a slab of normal conductor sandwiched between two superconductors (SNS junction). If we neglect the 
spatial dependence of the order parameter in superconductors, we don't need to solve the self-consistency equations 
(13110(1 . Therefore the problem reduces to a single-particle one and is most naturally solved in Bogoliubov-de Gennes 
language (Eq. becomes a Schrodinger equation for a two-component wave function) . 

Solutions of this equation with the boundary condition 112|) are standing waves. Obviously in the act of Andreev 
reflection a charge of ±2e is transferred to the superconductor, therefore every standing wave (Andreev level) carries 
supercurrent. 

Quasiclassically, in order to find Andreev levels in a normal layer of thickness L, sandwiched between "left" and 
"right" superconductors, with phase difference x, we write the Bohr-Sommerfeld quantization condition, 

I p(E)dq±x + Si(E)+6 r (E) = 2nn. (13) 

Here the kinematic phase gain of the quasiparticle along the closed trajectory, § p(E)dq — p e (E)dq + f r ph{E)dq = 
Ji (Pe(E) — Ph{E))dq, takes into account electron-hole (or vice versa) conversion. 

The positions of levels, and therefore the supercurrent, depend on the phase difference x between the supercon- 
ducting banks, and we arrive at Josephson effect in SNS structure a 5 ! 6 ' 7 . Actually, the language of Andreev levels can 
be successfully used to describe the Josephson effect in general (for a review seem). 



B. 7r-junctions and time-reversal symmetry breaking 

The crucial experiments(^i^*i^; see also review^) which confirmed d-wave pairing symmetry in high-T c cuprates 
were directed at catching the intrinsic phase 7r-shift. The general idea of the experiment follows from the fluxoid 
quantization condition in a superconductor: if a superconducting contour C is penetrated by the magnetic flux $, 
thenii 

$ / he 
27T — + & ds ■ = 27m, $ = — (14) 
$0 Jc 2e 

(in CGS units; $ ~ 2 • lCT 15 Wb in SI). 

In the case of a massive superconducting ring the contour can be chosen well inside the superconductor, where there 
is no current and therefore no superconducting phase gradient. Then the magnetic flux is quantized in units of $o- 

If there is a Josephson junction in the ring, the phase change will concentrate there, yielding 

2^ + X = 27m. (15) 

Here x is the phase difference across the junction. 

The equilibrium value of \ is determined by the interplay between Josephson and magnetic energy of the system. 
The Josephson energy is related to the Josephson current via 

KX) = 2e^|M, (16) 

and in the simplest case of a tunneling junction, I(x) — IcSinx, and E(\) — — (I c /2e) cosx. For a conventional 
Josephson junction E{\) is at a minimum when x — 0; therefore in the absence of an external field the total energy 
of the ring, 



U(x) = E( X )+E magn .^) = -^ C o SX +[^) gj, (17) 

has a single global minimum at \ = for any value of 2irLI c /$o (here L is the self-inductance of the ring). 

Not so if the pairing symmetry is <i-wave: if we choose the configuration of the loop in such a way, that the 
opposite sign lobes contact across the junction (so called 7r-junction), the current-phase relation switches to I(x) = 
7 c sin(x + 7r) = — J c sinx, and E(x) = +(I c /2e) cos X - Therefore if 2irLI c /$o > 1, the system has two degenerate 
minima; if this ratio is so big that the magnetic energy can be neglected compared to the Josephson energy, the 
equilibrium phase difference is xo = an d we obtain from (|15[) . that 

$=fn+^W (18) 
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Thus in equilibrium there is a spontaneous flux $o/2 in the ring. Its direction (up or down) allows us to distinguish 
the two ground states of the system, where the time-reversal symmetry is thus broken. 

Such behaviour was indeed observed in an experiment: a tri-crystal ring of YBCO generated a half-flux quantum, 
which was detected by SQUID microscopy. 

A natural question is whether 7r-junctions are the only possibility provided by rf-wave symmetry, and whether only 
$o/2 fluxes can be spontaneously generated. The answer is no: in principle, any equilibrium phase difference can be 
realized in a d-wave junction, and time-reversal symmetry breaking can be accompanied by generation of an arbitrary 
magnetic flux, or none at all. 



C. Josephson effect and T-breaking in SND and DND junctions 

Let us return to an SNS junction (assuming a rectangular normal part, L x W). Each quasiclassical trajectory 
connecting two superconductors acts as a conduit of supercurrent between them, an "Andreev tube" of diameter 
~ Af, which carries supercurrent, determined by the phase difference between its ends. If we neglect the normal 
scattering at NS interfaces, there is no mixing between different trajectories, since Andreev reflections simply reverse 
the velocity, sending the reflected particle along the same path. It can be shown, that if L ^> £0, every such trajectory 
passing through the point r contributes a partial supercurrent density^ 

Jl ' VF ' ~ X F W ^ [ ' l T smh(n£(v F )/l T f ' { J > 

In Eq.|JT!|J|, £(vp) is the length of the trajectory, and It — vp /2irksT is the so called normal metal coherence 
length. (We also included effects of weak impurity scattering with scattering time Tj.) The physical meaning of It 
is that in a clean normal metal an electron and Andrecv-rcflcctcd hole (or vice versa) with energy ksT maintain 
phase coherence across distance It simply because they travel along the same trajectory. Indeed, the momentum of 
an electron (hole) with energy ksT <C Ep is p e _h ~ Pf ± (dEPF)kBT = pp ± ksT/vp. At a distance I from the point 
of Andreev reflection (measured along the trajectory) they would gain phase difference (p e — Ph)l — 2ksTl/vF- If 
the phase difference is of order tt, the coherence is effectively lost, so we get /^ballistic ~ vp/ksT. The factor of 2tt 
appears in accurate treatment, like in Ea. (|19|l . 

As an aside, in case of very strong scattering in the normal part of the system, when the motion of electrons/holes 
is diffusive with diffusion coefficient D, we can still use the same argument. Now the observable length scale is 
given by the displacement of quasiparticle, I 2 = Dt, while the phase difference between the electron and the hole 
is gained along the crooked path of length I' — Vpt — vpl 2 /D they take. So the condition 2ksTV / 'vp ~ tt yields 

^T,diffusivc ~ \/ D/ksT. 

Such coherence in a normal metal is a purely kinematic effect, since there is no interaction in the normal (non- 
magnetic) metal, which would either support or suppress superconductivity. ("Normal metal is neutral with respect 
to superconducting correlations. "-C. W.J. Beenakker.) Nevertheless its effects are quite real, e.g. the very possibility 
of coherent supercurrent flow (Josephson effect) in SNS junctions. 

Let us apply the approach of 119|) to calculation of Josephson current in a clean SNS contact. In every point of 
NS boundary (e.g. x = L), we integrate l|19|> over directions of vp (such that vp x > 0). As it should be expected, 
the result reproduces Ishii's sawtooth^: in the limit of zero temperature and no scattering inside the normal layer 
(X = <h ~ fc), 

I(X) = [dy £ j x (y,v F ) oc 1 £)(-i)«+i***) = * W < n (20) 

/ * — ' TT * — ' U TT 

(periodically extended). 

Consider now an SND junction, that is, an SNS contact, where one of the superconductors has d-wave symmetry^ifi. 
Now, when integrating over the directions of vp we encounter two kinds of trajectories: "zero"- and "^"-trajectories, 
which link the conventional superconductor with the lobes of the <i-wave order parameter of opposite sign (intrinsic 
phase 7r-shift). The Josephson current is therefore a sum of two contributions like l|20ll . one of which has the phase 
argument x+7r. The relative weight of these contributions depends on the orientation of the SD boundary with respect 
to the crystal axes of the cuprate (to which the order parameter is nailed). In the case of a 45° orientation, when the 
boundary is along the (110) plane and is therefore normal to the nodal direction of the d-wave order parameter, the 
two groups contribute equally, yielding (cos 6* = vp x /vp) 

^ zeroeve s ^ cQg & ^ oo _ 2nL/l . cosf) L sin{n X ) + shl(n(x + 7r)) 

{X) ^ L tt ' ItcosO sinh(nL/ l T cos 6) ' 1 ' 

VFu:>0 n=l ' 
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FIG. 1: Josephson and spontaneous currents in SND and DND junctions. Arrows a and b indicate Andreev zero- and 7r-levels. 
In equilibrium their contributions to the Josephson current (normal to the NS boundary) cancel. By switching current directions 
(dashed arrows) we come to the other ground state, with the opposite direction of the spontaneous current, js. 




FIG. 2: Josephson and spontaneous currents in SND and DND junctions, a) Contributions of zero- and 7r-levels to the Josephson 
current in a 0-45° DD(SD) junction. b)Resulting current-phase dependence, c) Same in more general case, d) Spontaneous 
current in a 0-45° junction. 

All odd harmonics cancel, and we obtain a 7r-periodic sawtooth: at T = 0, li — > oo 

I(X)* X ~1 /2 , 0<X<k, (22) 

IT J A 

(periodically extended). (This is simply a sum of two identical, 27r-periodic sawtooth functions, shifted by 7r, see FigEl) 
There are two stable equilibrium phase differences across the junction: \ — ±7r/2. This means, that the time-reversal 
symmetry is broken. (In addition, the frequency of ac Josephson effect in the system will be doublediSiSl) . 

The SND junction may seem simply a "superposition" of two SNS junctions, with phase differences shifted by w, 
but the situation is more interesting. In equilibrium, spontaneous currents flow in the normal layer, parallel to the 
boundary^. This is clear from our "Andreev tube" treatment. Take, for example, two trajectories, with 9 = +a 
and 9 = —a, carrying equal partial currents. Their contribution to the Josephson current (normal to the boundary) 
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FIG. 3: Josephson and spontaneous currents in DND junctions—*—, a) Arbitrary case, b) Asymmetric (0-45°) junction, c) 
Symmetric (22.5°-22.5°) junction. Directions corresponding to zero-levels are shaded. In the case c) the contributions of zero- 
and 7r-levels to the current in y-direction are identically zero; therefore the spontaneous current is absent, but the time-reversal 
symmetry is still broken, d) Orientation of the d-wave order parameter with respect to the crystalline axes. 

is zero, since their currents' projections on this direction cancel each other. On the contrary, the projections on the 
direction, parallel to the boundary, add. The phase dependence of this current, I s (x)> can be obtained in the same 
way as for the Josephson current. In the limit of zero temperature and clean normal layer it reduces to a difference 
of two sawtooth functions, yielding 

I s (x) oc sgnx, |xl < T, (23) 

again periodically extended. States with 7r/2 and — 7r/2 thus carry spontaneous currents in opposite directions. If 
an SND junction is closed on itself (annular geometry), these currents translate into spontaneous magnetic fluxes, 
normal to the plane of the system^. 

Note that all of the above is only possible if higher harmonics of current-phase dependence are not all zero, since 
all the odd harmonics (including the standard Josephson term, sin \) cancel. (The expected " total depairing" in the 
(110) plane.) 

In the case of arbitrary orientation of the ND plane , the weights of zero- and 7r-sawtooth functions will be different 
(FigUp). Then the current-phase dependence regains 2-7r-periodicity. At zero temperature and in the absence of 
scattering, the T-breaking is generally still present. The equilibrium phase difference is no longer ±7r/2 and depends 
on the geometry. (The exception is when the boundary is along (100) or (010). Then the contribution of only one 
group survives, and we have either standard, or 7r-junction. Of course, it is impossible to tell, whether a single junction 
is a 7r-junction - it is necessary to look at the contour, in which it is included 12 .) Finite temperature and scattering 
suppress higher harmonics first, and therefore the transition to non- T-breaking state is possible^. 

Essentially the same interplay of zero- and 7r-levels takes place in DND junctions. There we have a somewhat richer 
picture. For example, the time-reversal symmetry can be broken without producing spontaneous currents (Figyp). 

One of the reliable methods of fabrication of Josephson devices in high-T c cuprates is based on forming grain 
boundary junctions (see e.g^). The order parameter is suppressed within ~ £o around the boundary, and this 
region can be considered "normal". Therefore the SND/DND model applies, but only qualitatively, since the above 
equations are derived in the limit L » (o. A more accurate approach, based on Eilenberger equations ©, confirms 
the qualitative similarity between DND and DD junctions^. 
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(a) (b) 




<J>/7T Cj>/7T 

FIG. 4: Josephson current (a) and spontaneous current (b) versus the phase difference in a clean DD grain boundary junction 
calculated in the non-self-consistent approximation. Temperature is T — 0.1T C . The mismatch angles are \i — and \ r = 45° 
(1), 40° (2), 34° (3), 22.5° (4)22. 

Integrating the current-phase dependences of FigEJ we obtain the Josephson energy of the junction, Ea. (|16|l . 
which has two minima, corresponding to its degenerate ground states. This bistability plays the crucial role in qubit 
applications of high-T c superconductors. This prediction was confirmed in grain boundary YBCO junctions^ (Fig|B). 

So far we did not take into account normal scattering on NS interfaces, and therefore missed an important point. 
Let us first consider a boundary of a d-wave superconductor with a vacuum or an insulator. The order parameter near 
the surface is suppressed, and a qualitative understanding of the situation can be obtained by "inserting" a normal 
layer of thickness ~ £o between the insulator and the bulk superconductor, similar to the DND model^t Energy levels 
in such a layer can be found from l|13|) for every quasiparticle trajectory (assuming specular reflection), FigJJp,: 

I p(E)dq + 26(E) + tts = 2nn. (24) 

Here s = or s — 1 depending on whether the trajectory connects the lobes of the same or opposite sign. Note that 
§ p(0)dq = and 26(0) — it. Therefore solution E = exists if and only if s = 1. This is the zero-energy, or midgap, 
state (ZES, MGS), which obviously cannot exist in conventional superconductors, where s always equalis zero (see 
review2£). 

Now consider a DD junction with finite transparency, T> < 1 (instead of the case of ideal transparency, T> = 1, 
which we dealt with earlier). We can use the DND model, inserting in the middle of the normal layer an infinitely 
thin barrier with transparency T>(6), dependent on the incidence angle of the quasiparticle trajectory^. 

Bohr-Sommerfeld quantization conditions for quasiparticle trajectories shown in Fig[7|3 yield the energy levels of 
the bound states in the junction. In the limit of low transparency, T> <C 1, the critical current is much larger (0(\fT>), 
not 0(T>)), if the orientations of the d-w&ve order parameter allow formation of zero energy states on both sides of 
the barrier. 

Even when transparency is not small, the presence of ZES may be qualitatively important. For example, the 
junction can become a 7r-junction at low enough temperature (see review^); you can see such behaviour in FigO), 
where below 11 K the potential wells are centered around it rather than zero. 

The appearance of spontaneous currents in SD and DD junctions can be considered as due to a time-reversal 
symmetry breaking order parameter with s + e lXo d or d + e lXo d' symmetry, which is formed in the junction area due 
to the proximity effect (here \o ^ 0, tt is the equilibrium phase difference across such junction). In certain conditions, 
such combinations could appear near the surface of a <i-wave superconductor, leading to spontaneous currents and 
magnetic moments. So far there is no conclusive evidence for such currents (see2^ and references therein). 
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5% 

FIG. 5: Equilibrium phase difference (f>o in a clean grain boundary junction as a function of the mismatch angle: \i = 0; 
X>- = (Fig- EJ- Parameter t = T/T c is the dimensionless temperature. Triangles and circles correspond to self-consistent and 
non-self-consistent solutions of Eilenberger equations, respectively. The 0-45° junction is a 7r/2-junction at any temperature, 
since by symmetry all odd components of the Josephson current are cancelled. Compare this to the behaviour observed in 
YBCO junctions (Fig®. 

II. MESOSOPIC d-WAVE QUBITS 
A. Flux qubits with conventional superconductors 

Qubits are the basic building blocks of future quantum computers. Essentially they are two-state quantum systems 
which can be put in an arbitrary superposition of states (" initialized" ) , coupled to each other, undergo desired quantum 
evolution and measured ("read out") before losing quantum coherence. Here we concentrate on superconducting phase 
qubits. 

The simplest example of such a qubit is an RF SQUID, that is, a loop with a single Josephson junction (like in 
Section ll B|l . We saw, that in certain circumstances the system has two degenerate minima, corresponding to a flux 
of ±$ /2 through the loop. 

The Hamiltonian of such a system is 

H=Uj(x) + U c (Q), (25) 

where Uq is the Coulomb energy of charge Q on the junction (which has some finite capacitance C). The charge 
operator Q can be expressed in terms of the phase difference across the junction as Q = —id x (see e.gAii). Due to the 
presence of the Coulomb term, phase is no longer a "good" quantum variable, and the eigenstates of the Hamiltonian 
(12511 become linear combinations of "up" and "down", or "left" and "right", states (with spontaneous flux ±<I>o/2). 
In other words, the qubit can now tunnel between the wells of the Josephson potential, corresponding to certain 
phases (FigEJ). (This description is appropriate as long as the Coulomb energy does not exceed the Josephson energy, 
otherwise the natural starting point would be the states with definite charge on the junction; we do not consider such 
systems (charge qubits) here, but they were successfully implemented experimentally 2 ^.) 

Coherent quantum tunneling was indeed observed in an RF SQUID 29 . Simultaneously, this effect was obtained in a 
different system, consisting of a small inductance loop with three Josephson junctions (the so called persistent current 
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FIG. 6: Bistability in a grain boundary YBCO junction—. Free energy, restored from Josephson current-phase dependence, is 
plotted vs. phase for (a) 0-45° and (b) symmetric (22.5°-22.5°) junction. Curves (top to bottom) correspond to temperatures 
(a) 30, 20, 15, 10, and 4.2 K; (b) 20, 15, 11, 10, 5, and 1.6 K, respectively. 




a) b) 



FIG. 7: Surface bound state in a d-wave superconductor (a) and model of a DD junction with finite transparency (b) (aftei— ) 



qubit^). 

The advantage of the latter design is as follows. As we have seen in Section ITbI the degenerate states appear in the 
RF SQUID only if the self-inductance of the loop is large enough, and they carry comparatively large spontaneous 
fluxes, ±$o/2. Therefore they will couple to the external degrees of freedom and reduce the time Td during which 
the system maintain its quantum coherence. Even more important is the fact that the resulting potential barrier is 
comparatively high, which may prohibit the tunneling we are after. In an experiment*^ the coherent tunneling was 
indeed observed not between the lowest, but between the excited states in the wells. 

In case of the three-junction loop of negligible self-inductance, the fluxoid quantization condition (|15|l leaves two 
independent Josephson phase differences in the circuit: 

Xi+X2+X2 + 2^f^ = 2 7 m. (26) 

The Josephson energy, 



Uj = -Eji cos(xi) - E J2 cos(x 2 ) - E J3 cos(x 3 ) (27) 
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FIG. 8: A tricrystal high-T c qubilA The ring D formed of high-T c film contains one 7r-junction and therefore supports a 
spontaneous flux ±<J?o/3"™- The ring S of conventional superconductor screens this flux from the environment. 

of the system forms a 2D potential profile (e.g. as a function of Xi, which depends on the external flux $ external = 
/ x $ as a parameter. 

If / = 0.5, and the -Ej's are comparable, this potential has two degenerate minima; unlike the case of the RF 
SQUID, the potential barrier can be small, as are the spontaneous fluxes corresponding to the "left" and "right" 
states. 

B. Rationale and proposed designs for qubits with d-wave superconductors 

One of the main problems with the designs of the previous section is the necessity to artificially break the T- 
symmetry of the system by putting a half flux quantum through it. Estimates show^i that the required relative 
accuracy is 10 -5 — 10 -6 . The micron-size qubits must be positioned close enough to each other to make possible their 
coupling; the dispersion of qubit parameters means that applied fields must be locally calibrated; this is a formidable 
task given such sources of field fluctuations as fields generated by persistent currents in qubits themselves, which 
depend on the state of the qubit; field creep in the shielding; captured fluxes and magnetic impurities. Moreover, the 
circuitry which produces and tunes the bias fields is an additional source of decoherence in the system. 

These problems are avoided if the qubit is intrinsically bistable. The most straightforward way to achieve this is 
to substitute the external flux by a static phase shifter, a Josephson junction with unconventional superconductors 
with nonzero equilibrium phase shift xo- For example, three-junction (persistent current, Mooij) qubits would require 
an extra 7r-junction (xo = 7r) 31 - The only difference compared to the case of an external magnetic field bias is in 
the prevalent decoherence sources: instead of noise from field-generating circuits we will have to take into account 
intrinsic decoherence from nodal quasiparticlcs and interface bound states (see below). 

Another suggestion^ is based on the same tricrystal high-T c ring geometry as in exp erimentsi&iii (Fi^SJ). The 
spontaneous flux ±$0 generated by such a structure labels the qubit states |0), |1). Tunnelling between them, 
necessary for quantum operations, is made possible by applying a magnetic field in the plane of the system. Indeed, 
then the states |0), |1) are no longer the eigenstates of the Hamiltonian. 

The tricrystal ring D is surrounded by an s-wave superconductor ring S, aimed at screening the spontaneous flux 
from the environment (including other qubits). Indeed, due to the fluxoid quantization condition (|14f> in the ring S, 
the total flux through it must be an integer multiple of $0, and therefore the states of the rings D and S become 
entangled (e.g. a\0) D |l) s + p\l) D |0) s ). 

This entanglement puts forward an interesting problem. In order to perform two-qubit operations, as well as 
initialization and readout, it is necessary to make the qubits "visible" to each other and the outside world. To do 
this, it is suggested ir>S to locally destroy superconductivity in the screening ring by using a superconducting field 
effect transistor (SUFET) (not shown) (i.e. by applying a gate voltage to the part of the screening ring). Will this 
transition collapse the wave function of the qubit, destroying quantum coherence between |0)d and |1)d? 

Now let us consider the case when the bistable d-wave system is employed dynamically, that is, when its phase 
is allowed to tunnel between the degenerate values. In a so called "quiet" qubili^ an SDS' junction (effectively two 
SD junctions in the (HO) direction) put in a small-inductance SQUID loop in parallel with a conventional Josephson 
junction and a large capacitor, Fig|5] One of the SD junctions plays the role of a 7r/2-phase shifter. The other 
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FIG. 9: "Quiet" SDS' qubit (aftej--). Switches c and s connect the SDS' structure to a conventional Josephson junction, S, 
and a large capacitor, C ox t. 

junction's capacitance C is small enough to make possible tunneling between the tt/2 and — 7r/2-states due to the 
Coulomb energy Q 2 jlC. Two consecutive SD junctions are effectively a single junction with equilibrium phases 
and 7r (which are chosen as working states of the qubit, |0) and |1)). The proposed control mechanisms are based 
on switches c, s. Switch c connects the small S'D junction to a large capacitor, thus suppressing the tunneling. 
Connecting s for the duration At creates an energy difference AE between |0) and |1), because in the latter case we 
have a frustrated SQUID with 0- and n- junctions, which generates $0/2 spontaneous flux. This is a generalization of 
applying a a z operation to the qubit. Finally, if switch c is open, the phase of the small junction can tunnel between 
and 7r. Entanglement between different qubits is realized by connecting them through another Josephson junction 
in a bigger SQUID loop. (The switches in question are in no way trivial, since they must operate without destroying 
quantum coherence of the system. One possible solution is to use a frustrated dc-SQUID 33 , that is, insert in the wire a 
small inductance loop with two equivalent Josephson junctions in parallel. The total supercurrent through the switch, 
/ = Ii(x%) +I2(X2), g° es to zero if the phase difference (tunable by external magnetic flux) \i ~ X2 — n - Modifications 
of this scheme are discussed m£i. Another possibility is to use superconducting single-electron transistors (SSETs, 
"parity keys")Si.) 

Due to the absence of currents through the loop during tunneling between |0) and |1) the authors called it "quiet", 
though as we have seen small currents and fluxes are still generated near SD boundaries. 

Another design based on the same bistability2iM only requires one SD or DD boundary fFig ll(J|) . Here a small 
island contacts a massive superconductor, and the angle between the orientation of the <i-wave order parameter and 
the direction of the boundary can be arbitrary (as long as it is compatible with bistability 21-22 , Fig|5J). The advantage 
of such a design is, that the potential barrier can be to a certain extent controlled and suppressed; moreover, in 
general the two "working" minima (— </>o, 4>o, the phase of the bulk superconductor across the boundary is zero) will 
be separated from each other by a smaller barrier, than from the equivalent states differing by 2irn. This allows us to 
disregard the "leakage" of the qubit state from the working space spanned by (|0),|1)), which cannot be done in the 
"quiet" design with the exact 7r-periodicity of the potential profile. 

Qubit operations in this system are realized by connecting qubits to each other and to the ground electrode (normal 
or superconducting) through SSETs (or other kind of switches). When isolated, a qubit undergoes natural evolution 
between |0) (state with phase — <fio) and |1) (phase <fio), which realizes the a x operation. The a z operation (that is, 
adding a controllable phase shift to one of the states with respect to the other) can be realized by e.g. connecting 
the island through a SSET to the massive superconductor ("bus"), the phase of which The same operation 

repeated periodically can be used to block unwanted tunneling (so called "bang-bang" technique)^: physically, if we 
keep shifting the levels in the right and left well with respect to each other, the tunneling becomes suppressed, since 
they are practically never in resonance. 

A better design was suggested inS ("silent qubit"). Here two small bistable <i-wave grain boundary junctions with 
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FIG. 10: A qubit based on an SD (a) or grain boundary DD (b) Josephson junction—*—. Parity keys (superconducting 
single electron transistors) are used to connect the nearest neighbours. The critical current through the parity key (a small 
superconducting island) can be tuned by changing the gate voltage, V^— . M are magnetic force microscope tips, intended to 
read out spontaneous fluxes. An alternative readout scheme relies on direct detection of the phase on the island, which affects 
the critical current between the bus and the ground—. 



a small superconducting island between them are set in a SQUID loop (Fig. Ill|) . As usual, "small" means that the 
total capacitance of the system allows phase tunneling: the Coulomb energy term in 1|25|) is not negligible. 

In the limit of negligible self-inductance of the loop, the quantization condition l|14l) fixes the sum of the phases to 
the external flux, \i + X2 = 4> = 27r$/<I>o- This leaves only one independent phase combination, the superconducting 
phase of the island, = (xi ~ X2)/2- 

Keeping for simplicity only the first two harmonics in the current-phase relation of the junctions, 

h = I> smxi - I- sin 2\ = 
= I 0i [sin Xi sin 7* - sin 2\i cos 7,-] , (28) 

we find for the Josephson potential of the qubit the expression (Fig. I12f> 

Uj(6, 0) = -(/ i/2e) [f(<j>/2 + 0, 7l ) + v f(fc - 0, j 2 )] . (29) 

Here f(<p, 7) = cos^) sin(7) — (l/2)cos(2ip) cos(7), and r\ = I^/Iqi- Parameter 7 € [0, 7r/2] provides a convenient 
parametrization for the current-phase dependence in a DD junction. 

In the absence of the external flux (<j) = 0) the qubit potential Uj(0, 0) has two degenerate minima. Moreover, if the 
junctions only differ in the amplitude of critical current, but have the same 7, there is no spontaneous current in the 
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FIG. 11: A silent d-wave qubflA The structure is formed of a high-T c film around the grain boundary AA' and contains two 
grain boundary junction (J, J'). 




FIG. 12: The Josephson energy profile for the silent qubit^ as a function of the island phase, 8; 71 = 72 = 7r/4; 77 = 0.5. In 
the limit of negligible self-inductance the potential depends on the external magnetic flux, $ cx t = $0 x (<j>/2n) as on parameter 
(a). In the absence of external flux the potential has two degenerate minima; degeneracy is lifted by finite external flux, which 
also affects the height of the potential barrier(b). 

loop in either minimum, which means that the qubit is decoupled from the external magnetic fields (if we disregard 
spontaneous currents in the junctions themselves, which can be very small 21 *^^* 2 ^) . This justifies the moniker "silent" . 

Both the barrier and the bias between the wells can be controlled by the external flux. It is noteworthy that the 
corrections are of at least second order in <p, which drastically reduces the influence of fluctuations in the external 
circuitry^. The mechanism of noise reduction is similar to that of the " quantronium" qubili^, but the " sweet spot" 
(an extremal point on the energy surface) appears already on the classical level, and at zero external field. 

Finally, let us briefly mention two more proposals. 

A "no tunneling" design 40 combines the ideas of CBJJ (current-based Josephson junction) qubits 4142 and the 
intrinsic bistability of d-wave junctions. In a single bistable Josephson junction, the potential barrier between the 
degenerate levels corresponding to phase difference ±Xo is too high to allow tunneling. The transitions between the 
states are realized through Rabi transitions via an auxiliary energy level, situated above the top of the barrier. Rabi 
transitions are induced by applying an external high-frequency field. 

"Dot/antidot" proposals are based on the spontaneous flux (less than a flux quantum) generated in a high-T c island 
or around a hole in bulk high-T c due to the presence of a subdominant order paramete r 23 ! 36 (see the end of Section 
II C|) . but so far lack experimental support. 
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FIG. 13: Critical currents of two nominally identical dc SQUIDs, based on YBCO grain boundary junctions, as a function of 
applied magnetic field—. Dots: experiment. Thin line: theoretical fit. Fitting parameters: (a) J-f 7 ' = 9 fiA; J-f 77 ' — 3.7 /j,A; 
4 J) = 0.3 ix A; = 22.7 fxA; (b) = 7.8 /iA; l[ u) = 5.3 fiA; = 3 fiA; I { 2 n) = 4.3 fiA. 

C. Fabrication and experiment 

Due to fabrication difficulties, as well as expected problems with decoherence from e.g. nodal quasiparticles (see 
Section |11D|I . experimental research on d-wave qubits is not as far advanced as on the devices with conventional 
superconductors. The experimental confirmation of quantum behaviour in these systems is still missing. Nevertheless 
several recent successes should be noted. 

Arrays of half-flux quanta were realized and manipulated (classically) in YBCO-Nb zigzag junctions— Each facet 
of the junction effectively constitutes a 7r-ring, which supports spontaneous flux of ±$ - Interaction between the 
fluxes is mainly due to the superconducting connection, which leads to robust antiferromagnetic ordering. In its 
absence, a weaker, magnetic interaction establishes ferromagnetic flux ordering. The authors consider the possibility 
of using their structures in qubit design. 

Good quality submicron grain boundary YBCO junctions were fabricated and bistable energy vs. phase dependence 
was demonstrated 26-43 . Dc SQUIDs YBCO (15x15 /im 2 square loops with nominally 2 fim wide grain boundary 
junctions) were fabricated and tested, and their classical behaviour is very well described by the existing theorjaa (Fig. 
I13p . Like in Q28p. only two harmonics of current-phase dependence were considered. From fitting the experimental 
data, we had to conclude that the junctions in the same SQUID have not only different critical current amplitudes, but 
different ratios of first to second harmonic (71 7^ 72), probably due to the variation in the grain boundary properties 
over a distance of ~ 15 /im. This is one reason for putting the junctions in the silent qubit, Fig. 1111 close to each 
other, as it was done when fabricating its prototype— 

D. Decoherence in d-wave qubits 

Decoherence is the major concern for any qubit realization, especially for solid state qubits, due to abundance 
of low-energy degrees of freedom. In superconductors, this problem is mitigated by the exclusion of quasiparticle 
excitations due to the superconducting gap. This explains also why the very fact of existence of gapless excitations 
in high-T c superconductors long served as a deterrent against serious search for macroscopic quantum coherence in 
these systems. An additional source of trouble may be zero-energy states in DD junctions. 

Nevertheless, recent theoretical analysis of DD junctions 4 ^ 6 ^^, all using quasiclassical Eilenberger equations, shows 
that the detrimental role of nodal quasiparticles and ZES could be exaggerated. 

Before turning to these results, let us first do a simple estimate of dissipation due to nodal quasiparticles in bulk 
d-wave superconductors. 

Consider, for example, a three-junction ("Mooij") qubit with d-wave phase shifters. The |0) and |1) states sup- 
port, respectively, clockwise and counterclockwise persistent currents around the loop, with superfluid velocity v s . 
Tunnelling between these states leads to nonzero average (v 2 ) in the bulk of the superconducting loop. 

Time-dependent superfluid velocity produces a local electric field 



(30) 
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and quasiparticle current j qp — crE. The resulting average energy dissipation rate per unit volume is 

£ = aE 2 w mT qp (n{v s )v 2 s ). (31) 

Here T qp is the quasiparticle lifetime, and 

/>oc 

n(v s ) = / deN(e)[n,F(e- pfv s ) + riF(e+pFV s )} (32) 
Jo 

is the effective quasiparticle density. The angle-averaged density of states inside the d-wave gap is^ 

#( e ) „#(()) Jf_, (33) 

where pL — d ^g ^ > ancl is the maximal value of the superconducting order parameter. Substituting l|33|) in 13211 , 
we obtain 

n( Vs ) „ iV(0)-|- (-T 2 ) (Li 2 (-e- 2 ^) + Li 2 (-e 2 ^)) , (34) 

where Li2(z) = J eft ln(i-t) ^ g ^logarithm. Expanding for small pfv s <C T, and taking into account that 
Li 2 (-1) = -tt 2 /12! Li 2 (-1) = ln2, and Li 2 '(-1) = -1/2 + hi 2, we obtain 

The two terms in parentheses correspond to thermal activation of quasiparticles and their " Cherenkov" generation by 
current-carrying state. Note that finite quasiparticle density does not lead by itself to any dissipation. 
In the opposite limit (T <C pfv s ) only the "Cherenkov" contribution remains, 

n { v s) ~ 7— (PFVs) (36) 

(since Li 2 (z) ~ — (l/2)(lnz) 2 for large negative z, and ~ z for small z). 

The energy dissipation rate provides the upper limit r e for the decoherence time (since dissipation is sufficient, but 
not necessary condition for decoherence). Denoting by I c the amplitude of the persistent current in the loop, by L 
the inductance of the loop, and by £1 the effective volume of d-wave superconductor, where it flows, we can write 

i 2 s n 2 m r qp N(0)n[^f(v 2 )+p 2 F (v 2 v 2 )) 
Te ~ LP C ~ mAoLJ* ■ (37) 

Note that the thermal contribution to r" 1 is independent on the absolute value of the supercurrent in the loop (oc v s ), 
while the other term scales as I 2 . Both contributions are proportional to fl and (via v s ) to u>t, the characteristic 
frequency of current oscillations (i.e. tunneling rate between clockwise and counterclockwise current states). 

It follows from the above analysis that the intrinsic decoherence in a d-wave superconductor due to nodal quasipar- 
ticles can be minimized by decreasing the amplitude of the supercurrent through it, and the volume of the material 
where time- dependent supercurrents flow. From this point of view the designs with phase shifters, as well as the 
Newns-Tsuei and "no tunneling" designs are at a disadvantage (the latter, because the microwave field necessary to 
produce Rabi transitions will affect the whole sample). 

Now let us estimate dissipation in a DD junction. First, followingi2i2i, consider a DND model with ideally transmis- 
sive ND boundaries. Due to tunneling, the phase will fluctuate, creating a finite voltage on the junction, V — (l/2e)x, 
and normal current /„ = GV. The corresponding dissipative function and decay decrement are 

1* G± 2 fl^ 2 



T=-E = -GV = — ! (38) 

2 2 2 \2e) y 1 

2 dT G AN±En , x 

- Q . (39) 



M Q x d X 4e 2 M Q 
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Here Eq = e 2 /2C, Mq = C/16e 2 = 1/32Eq, N± are the Coulomb energy, effective "mass" and number of quantum 
channels in the junction respectively. The latter is related to the critical Josephson current Iq and spacing between 
Andreev levels in the normal part of the system e = vp/2L via 

J = N ± ee. (40) 

We require, that -f/u>o <C 1, where luq = y/32N±EQi/ir is the frequency of small phase oscillations near a local 
minimum. This means, 

N ± « JL (41) 

The above condition allows a straightforward physical interpretation. In the absence of thermal excitations, the only 
dissipation mechanism in the normal part of the system is through the transitions between Andreev levels, induced 
by fluctuation voltage. These transitions become possible, if e < 2eV ~ \fjc* ~ w , which brings us back to Q4ip. 
Another interpretation of this criterion arises if we rewrite it as lo^ 1 3> [vp / L)~ 124 . On the right-hand side we see 
time for a quasiparticle to traverse the normal part of the junction. If it exceeds the period of phase oscillations 
(on the left-hand side), Andreev levels simply don't have time to form. Since they provide the only mechanism for 
coherent transport through the system, the latter is impossible, unless our "no dissipation" criterion holds. 

For the thickness of the normal layer L ~ 1000 A and vf ~ 10 7 cm/s this criterion limits usq < 10 12 s _1 , which 
is a comfortable two orders of magnitude above the usually obtained tunneling splitting in such qubits (~ 1 GHz) 
and can be accommodated in the above designs. Nevertheless, while presenting a useful qualitative picture, the DND 
model is not adequate for the task of extracting quantitative predictions. For example, the coherence length It in 
the normal metal can be very large, while in the high-T c compound it is short. Therefore the estimates for crucial 
parameters (like e) based on the assumption It 3> L can be wrong. Moreover, the assumption of ideal ND boundaries 
is not realistic. 

A calculation^, which used a model of a DD junction interacting with a bosonic thermal bath, gave an optimistic 
estimate for the quality of the tricrystal qubit Q > 10 s . 

The role of size quantization of quasiparticles in small DD and SND structures was suggested in2ii£. The importance 
of the effect is that it would exponentially suppress the quasiparticle density and therefore the dissipation below 
temperature of the quantization gap, estimated as 1 — 10 K. Recently this problem was investigated^ for a finite 
width DD junction. Contrary to the expectations, the size quantization as such turned out to be effectively absent 
on the scale exceeding £o (that is, practically irrelevant). On the other hand, the finite transverse size of the system 
imposed an effective band structure. Namely, due to the direction dependence of the order parameter, a quasiparticle 
travelling along the trajectory, bouncing between the sides of the junction, goes through a periodic ID off-diagonal 
potential. The influence of this band structure on dissipation in the system is not straightforward. Ironically, from 
the practical point of view this is a moot point, since the decoherence time from the quasiparticles in the junction, 
estimated inM^, already corresponds to the quality factor T v /r g ~ 10 6 , which exceeds by two orders of magnitude 
the theoretical threshold allowing to run a quantum computer indefinitely. 

The expression for the decoherence time obtained in 45 ' 4§ , 



4e 



8<hI(At/ey 



(42) 



where 5<p is the difference between equilibrium phases in degenerate minima of the junction (i.e. 5<p = 2\o in other 
notation), contains the expression for quasiparticle current in the junction at finite voltage At/e (where At is the 
tunneling rate between the minima). This agrees with our back-of-the-envelope analysis: phase tunneling leads to 
finite voltage in the system through the second Josephson relation and with finite voltage comes quasiparticle current 
and decoherence. The aforementioned quality factor is defined as Q — T v A t /2h, that is, we compare the decoherence 
time with the tunneling time. This is a usual optimistic estimate, since it will certainly take several tunneling cycles 
to perform a quantum gate operation. 

It turns out that a much bigger threat is posed by the contribution from zero energy bound states, which can be at 
least two orders of magnitude larger. We can see this qualitatively from l|42[l : a large density of quasiparticle states 
close to zero energy (i.e., on Fermi level) means that even small voltages create large quasiparticle currents, which sit 
in the denominator of the expression for t v . Fortunately, this contribution is suppressed in the case of ZES splitting, 
and such splitting is always present due to, e.g., finite equilibrium phase difference across the junction^. 

A similar picture follows from the analysis presented in^S,. A specific question addressed there is especially important: 
it is known that i?C-constant measured in DD junctions is consistently 1 ps over a wide range of junction sizes^, 
and it is tempting to accept this value as the dissipation rate in the system. It would be a death knell for any 
quantum computing application of high-T c structures, and nearly that for any hope to see there some quantum 
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FIG. 14: Phase (a) and temperature (b) dependence of the normal conductance of a DD junction— i. 

effects. Nevertheless, it is not quite that bad. Indeed, we saw that ZES play a major role in dissipation in a DD 
junction, but are sensitive to phase differences across Measurements of the RC constant are done in the resistive 
regime, when a finite voltage exists across the junction, so that the phase difference grows monotonously in time, 
forcing ZES to approach the Fermi surface repeatedly. Therefore trc reflects some averaged dissipation rate. On the 
other hand, in a free junction with not too high a tunneling rate phase differences obviously tend to oscillate around 
Xo or — xoi its equilibrium values, and do not spend much time near zero or tt; therefore ZES are usually shifted from 
the Fermi level, and their contribution to dissipation is suppressed. 

This qualitative picture is confirmed by a detailed calculation^ (Fig. 114(1 . The decoherence time is related to the 
phase-dependent conductance via 



1 , SE 

t,„ = : — -77— — tanh . 

v aF{ Xo ) 2 SE 2T 

Here a is the dissipation coefficient, 5E is interlevel spacing in the well, and 

G( X ) = 4e 2 a {d x F( X )f . 



(43) 



(44) 



For a realistic choice of parameters Ea. (|43|l gives a conservative estimate t v = 1 — 100 ns, and quality factor Q ~ 1—100. 
This is, of course, too little for quantum computing, but quite enough for observation of quantum tunneling and 
coherence in such junctions. 



III. CONCLUSIONS 



Since the d-wave character of superconducting pairing in high-T c cuprates was established, our understanding of 
high-T c Josephson structures and ability to fabricate them progressed significantly. Now it enables us to ask the 
question in the title of this paper: Can high-T c cuprates play a role in quantum computing?- and to tentatively 
answer: Yes. 

We have seen that submicron junctions of sufficiently high quality are now fabricated. Several designs, which take 
advantage of the intrinsic bistability of d-w&ve structures, were developed, fabricated, and tested in classical regime. 
If the decoherence time turns out to be large enough, their better scalability can be a decisive advantage for quantum 
computing applications. 

This is, of course, a big if. Still, theoretical estimates tend to be rather optimistic. Even though they widely 
differ, all of them predict t v long enough to encourage experimental search for quantum coherence. There are several 
compelling reasons to do that. First, it would be a spectacular result. Second, it would clarify why different models 
give different answers. Third, it could indeed lead to practical application of high-T c devices in quantum computing. 
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P.S. 

Since this paper was published (2003), macroscopic quantum tunneling was observed in YBCO*& and BiSCCO^l 
junctions. This made the "big if of the previous paragraph somewhat smaller. 
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